Simulations of the Nonlinear Helmholtz Equation: Arrest of Beam Collapse, Nonparaxial Solitons, and Counter-Propagating Beams

We solve the (2+1)D nonlinear Helmholtz equation (NLH) for input beams that collapse in the simpler NLS model. Thereby, we provide the first ever numerical evidence that nonparaxiality and backscattering can arrest the collapse. We also solve the (1+1)D NLH and show that solitons with radius of only half the wavelength can propagate over forty diffraction lengths with no distortions. In both cases we calculate the backscattered field, which has not been done previously. Finally, we compute the dynamics of counter-propagating solitons using the NLH model, which is more comprehensive than the previously used coupled NLS model.Comment: 6 pages, 6 figures, Lette

A High-Order Numerical Method for the Nonlinear Helmholtz Equation in Multidimensional Layered Media

We present a novel computational methodology for solving the scalar nonlinear Helmholtz equation (NLH) that governs the propagation of laser light in Kerr dielectrics. The methodology addresses two well-known challenges in nonlinear optics: Singular behavior of solutions when the scattering in the medium is assumed predominantly forward (paraxial regime), and the presence of discontinuities in the % linear and nonlinear optical properties of the medium. Specifically, we consider a slab of nonlinear material which may be grated in the direction of propagation and which is immersed in a linear medium as a whole. The key components of the methodology are a semi-compact high-order finite-difference scheme that maintains accuracy across the discontinuities and enables sub-wavelength resolution on large domains at a tolerable cost, a nonlocal two-way artificial boundary condition (ABC) that simultaneously facilitates the reflectionless propagation of the outgoing waves and forward propagation of the given incoming waves, and a nonlinear solver based on Newton's method. The proposed methodology combines and substantially extends the capabilities of our previous techniques built for 1Dand for multi-D. It facilitates a direct numerical study of nonparaxial propagation and goes well beyond the approaches in the literature based on the "augmented" paraxial models. In particular, it provides the first ever evidence that the singularity of the solution indeed disappears in the scalar NLH model that includes the nonparaxial effects. It also enables simulation of the wavelength-width spatial solitons, as well as of the counter-propagating solitons.Comment: 40 pages, 10 figure

Statistical Mechanics of Charged Particles in the Pressure of Magnetic Irregularities

Statistical mechanics of charged particles in presence of magnetic irregularitie

Ring-type singular solutions of the biharmonic nonlinear Schrodinger equation

We present new singular solutions of the biharmonic nonlinear Schrodinger equation in dimension d and nonlinearity exponent 2\sigma+1. These solutions collapse with the quasi self-similar ring profile, with ring width L(t) that vanishes at singularity, and radius proportional to L^\alpha, where \alpha=(4-\sigma)/(\sigma(d-1)). The blowup rate of these solutions is 1/(3+\alpha) for 4/d\le\sigma<4, and slightly faster than 1/4 for \sigma=4. These solutions are analogous to the ring-type solutions of the nonlinear Schrodinger equation.Comment: 21 pages, 13 figures, research articl

High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension

The nonlinear Helmholtz equation (NLH) models the propagation of electromagnetic waves in Kerr media, and describes a range of important phenomena in nonlinear optics and in other areas. In our previous work, we developed a fourth order method for its numerical solution that involved an iterative solver based on freezing the nonlinearity. The method enabled a direct simulation of nonlinear self-focusing in the nonparaxial regime, and a quantitative prediction of backscattering. However, our simulations showed that there is a threshold value for the magnitude of the nonlinearity, above which the iterations diverge. In this study, we numerically solve the one-dimensional NLH using a Newton-type nonlinear solver. Because the Kerr nonlinearity contains absolute values of the field, the NLH has to be recast as a system of two real equations in order to apply Newton's method. Our numerical simulations show that Newton's method converges rapidly and, in contradistinction with the iterations based on freezing the nonlinearity, enables computations for very high levels of nonlinearity. In addition, we introduce a novel compact finite-volume fourth order discretization for the NLH with material discontinuities.The one-dimensional results of the current paper create a foundation for the analysis of multi-dimensional problems in the future.Comment: 47 pages, 8 figure

Stable self similar blow up dynamics for slightly L^2 supercritical NLS equations

We consider the focusing nonlinear Schr\"odinger equations $i\partial_t u+\Delta u +u|u|^{p-1}=0$ in dimension $1\leq N\leq 5$ and for slightly $L^2$ supercritical nonlinearities p_c with $p_c=1+\frac{4}{N}$ and 0<\e\ll 1. We prove the existence and stability in the energy space $H^1$ of a self similar finite time blow up dynamics and provide a qualitative description of the singularity formation near the blow up tim

Singular solutions of the L^2-supercritical biharmonic Nonlinear Schrodinger equation

We use asymptotic analysis and numerical simulations to study peak-type singular solutions of the supercritical biharmonic NLS. These solutions have a quartic-root blowup rate, and collapse with a quasi self-similar universal profile, which is a zero-Hamiltonian solution of a fourth-order nonlinear eigenvalue problem